Two-mode variational Monte Carlo study of Quasiparticle excitations in cuprates 
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Recent measurements of quasiparticles in hole-doped cuprates reveal highly unusual features: 1) 
the doping-independent Fermi velocity, 2) two energy scales in the quasiparticle spectral function, 
and 3) a suppression of the low energy spectral weight near the zone center. The underlying mech- 
anism is under hot debate. We addressed these important issues by a novel two-mode variational 
Monte Carlo (VMC) study of the t- J model. We obtained results in agreement with the experiments 
but without invoking extrinsic effects. Besides, we resolved a long standing issue of the sum rule for 
quasiparticle spectral weights in VMC studies. The electron doped case was also discussed. 
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Understanding the quasiparticle excitations is a pre- 
clude to unravel the mechanism of high-Tc supercon- 
ductivity. Recent angle-resolved photo-emission spectra 
(ARPES) in cuprates reveal highly unusual features and 
lead to hot debates. For example, on the low energy 
quasiparticle band the nodal Fermi velocity does not 
seem to increase with increasing doping, and is called 
a "universal nodal Fermi velocity" . [1] This seems to 
challenge the concept of doped Mott insulators which 
would naively predict a vanishing Fermi velocity with 
decreasing doping. On the other hand, ARPES mea- 
surements up to higher energies [2-9] revealed that the 
quasi-particle dispersion along the nodal direction breaks 
up near the momentum (7r/4, 7r/4) at an energy around 
0.3 ~ 0.4 eV (or the low energy spectral weight near 
the F seems to disappear), and reappeared around leV. 
The two bands are connected in a waterfall fashion in the 
momentum- energy space. This high energy anomaly has 
intrigued many theoretical studies and disputes. [2, 10-17] 
A central concern is whether the anomaly is caused by 
extrinsic effects due to phonons, or is an intrinsic prop- 
erty of the strongly correlated electrons themselves. We 
summarize the main important issues that will concern 
us in this Letter: 1) the doping-independent Fermi ve- 
locity, 2) two energy scales in the quasiparticle spectral 
function, and 3) a suppression of the low energy spectral 
weight near the zone center. 

We address the above important issues within the com- 
monly accepted one-band t-J model without invoking 
other extrinsic effects. The theoretical machinery we 
use is the Gutzwiller projection variational Monte Carlo 
(VMC). It proves to give good energy and supports d- 
wave pairing symmetry. [18-22] As for quasiparticle ex- 
citations under concern, a projected mean-field excited 
state is commonly used in the literature. [23-25] For a 
given momentum k and spin ct, only one such quasipar- 
ticle can be constructed. We call such an approach as 
a single-mode approach (SMA). Unfortunately, the SMA 
gives only a single low energy band and predicts that 
the spectral weight is maximal at the center of the Bril- 
louine zone (the F point), in contrast to the two energy 
scales and suppression of low energy spectral weight re- 



vealed by ARPES. The spectral weight below the Fermi 
level captured by SMA is proved to be smaller than the 
exact value 1 — a;. [23, 26] It is therefore likely that the 
high energy feature is due to the missing spectral weight 
beyond the SMA. Moreover, the nodal Fermi velocity in- 
creases significantly with doping in such an approach. [25] 
In this Letter we design a two-mode approach (TMA) 
for each set of quasiparticle quantum number (k, a) . We 
show that the TMA satisfies the spectral sum rules for 
quasiparticle excitations both below and above the Fermi 
level, resolving a long standing problem SMA. [27] In this 
way spectral weights beyond the scope of the SMA can 
also be captured. Our numerical results are in agree- 
ment with the main experimental features listed above. 
We also demonstrate that in the electron doped case 
there should be no high energy anomaly below the Fermi 
level, in agreement with the finite temperature exact 
diagonalization. [14] 



The t-J model Hamiltonian is H = Ht + Hj, with Ht = 
-Efe),.*y(J^GC,Lc„.PG + h.c.) and Hj = JEfe)('5, • 
Sj - jPGntnjPG)-[28] Here Uj = ti,t2,h are hopping 
integral between the nearest, second- nearest and third- 
nearest neighbor sites i and j, Cia is the electron annihi- 
lation operator. Si is the electron spin, rii = mi + Uii 
with = c\^Cia, and Pq = 11^(1 - ni^mi) is the 
Gutzwiller projection operator that removes any dou- 
ble occupations. A widely used trial ground state with 
d-wave pairing is \^gs) = PgPn\^ dBCs) ,[^8] where 
Pn is the projection operator that fixes the number 
of electrons, J^dscs) = Ilfcl^fe + ^*:4T^-/ci)|0)- Here 

< = W + t)> - 5(1 - t)' ^'^ - v/f + ^' = 
— 2(cos kx + cos ky) — 4i^jj^ cos k^ cos ky — 2ty^^{cos 2kx + 

C0s2ky) l-Lvarj — ^varicOS kx COS ky) . t^^^,, i^^j., 

Ayar and /i^ar are variational parameters. In the follow- 
ing we concentrate on the quasiparticle excitations. Can- 
didates for such excitations are described by the quasi- 
hole, bare-hole, quasi-electron and bare-electron wave 
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functions as follows: 



given by 



\'^qh{k,a)) = 

|*6fc(fc,C7)) = 

|*ge(A;,cr)) = 



PcPN-lllj'^dBCs), 
C-k-al'^Gs), 

PgPn+hU^cIbcs), 



(1) 

(2) 

(3) 

(4) 



Here jka = UkCka- — crufcc^j,- is the Bogoliubov quasi- 
particle annihilation operator with momentum k and spin 
a. As pointed out by Ran et al, all these wave functions 
on for a specific real-space electron configuration can be 
written in the form of determinants. [29] Thus the over- 
laps between these states can be computed statistically. 
To ease further discussion we denote the normalized wave 
functions for the ground state, the quasi-hole, the bare- 
hole, the quasi-electron and the bare-electron states as, 
respectively, \GS), \QH), \BH), \QE) and \BE). The 
subscripts k and a arc left implicit. 

The ARPES experiment measures the one particle 
spectral function A{k,uj) in the occupied side (w < 0). 
According to the Lehnmann representation, 

A{k,Lo) =^[|H4 J0)|2<5(u; + u;o - u^n) 

n,a 

+ \{n\c-k,-a\0)\^S{w-wo+Un)]. (5) 

Here |n) denotes an excited eigenstate oi H — with 
the eigenvalue a;„. Wc define Z~ = \{n\cka\^)\'^ and 
■^n ~ K'^klo-I^)!^ ^ spectral weights of \n) in the oc- 
cupied and unoccupied sides, respectively. For a free par- 
ticle system, the state with the quantum numbers (k, a) 
is unique. This is no longer the case in an interacting sys- 
tem. We recall that \QH) is the only state used in SMA. 
As mentioned above, the spectral weights not captured 
by SMA are likely to appear at higher energies, and this 
motivates us to go beyond the SMA by enforcing the sum 
rule. 

In principle one should construct a complete set of ex- 
citations to satisfy the sum rule. This is possible by 
exact diagonalization but is limited by the lattice size. 
We therefore take a simpler route but still keep the 
sum rule. The idea is as follows. Since [cfco-j-Pc] 7^ 0, 
\QH) and \BH) arc neither identical nor orthogonal. 
We can reconstruct two orthogonal states \BH) and 
--{\QH) — k\BH)) for each set of quantum 



numbers (k, cr), where k = {BH\QH). The Hamiltonian 
matrix in this basis space is 



_/ {rj\H\n) {v\H\BH) \ 
{ii)2x2 - 1 i^BH\H\'n) {BH\H\BH) ) ' 



(6) 



which we diaganolize to get the eigenstates |1) and |2) as 

linear combinations of \BH) and |?7) (and thus of \BH) 
and \QH). The occupied spectral weight for (k, a) is now 



_ l(l|c-fc,- 



I* 



GSI 



^ |(2|c_fc,_.|vl'Gs)| 

(vl/fa^|(|l)(l| + |2)(2|)|vl/fe;,) 



(*GS|*GS) 



(*GS|*GS) 



k,—ai 



(7) 



where we used definition of the unnormalized bare-hole 
state. We also used the crucial fact that I'^bh) lies within 
the space spanned by |1) and |2), or equivalently by \BH) 
and 1 77). By summing over (k, a) we get a total occupied 
spectral weight J2ka ^ka/S = 1 — x, with S the number 
of lattice size, as required by the sum rule. [26] We call 
the above procedure as a two-mode approach (TMA) for 
the excitations in the occupied side (or below the Fermi 
level). On the unoccupied side, since Pg[c',Pg] ~ 
the bare-electron and quasi-electron states are identical 
within the Hilbert space of the model, and one can prove 
that the total spectral weight of \QE) already satisfies 
the sum rule J2ka Ka/^ = 2a;.[25, 26] 

We perform Monte Carlo calculations for the energy 
and spectral weight of the states |1) and |2) below the 
Fermi level, and of the state \QE) above the Fermi level. 
We also present the results under SMA for comparison. 
The calculation is done on a 10 x 10 lattice. Larger sizes 
up to 14 X 14 are also attempted which do not alter the 
conclusions we reach. In our calculation we use 200,000 
samples to reduce the statistical error, which we find to 
be less than 5 meV for quasiparticle energies. In or- 
der to get finer resolution in momentum space we use 
four combinations of periodic and anti-periodic boundary 
conditions. It turns out that the quasi-hole and quasi- 
electron bands are symmetric in energy, [25] and we define 
the central energy as the Fermi level. In Figs.l. the linear 
size of the symbols represents the momentum-dependent 
spectral weight, while the central position of the symbols 
represent the quasi-particle energy as a function of mo- 
mentum along the cuts (0, 0) (tt, tt) (tt, 0) (0, 0). 
The doping levels are x = 6% (a), x = 20% (b), and 
X = 30% (c). The fiUed-blue circles are for excitations on 
the unoccupied side, which are identical under both SMA 
and TMA. The open-green circles corresponds to the re- 
sult given by the SMA, which yields a single band below 
the Fermi level with the largest spectral weight at the F- 
point. This is a common result under SMA, [23, 24] but 
is inconsistent with the recent ARPES measurements. 
This should be contrast to the following TMA results. 
The fiUed-red and filled-pink circles are for excitations 
on the occupied side contributed by the lower energy 
eigenstate )1) and higher energy eigenstate ]2), respec- 
tively. We find that along the nodal direction the spectral 
weight of |1) state at the F point is greatly suppressed 
and gradually increases away from the F point in the 
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FIG. 1: (Color online) Energy dispersion and spectral weight 
map in the energy-momentum space (along high symmetry 
cuts). Only nearest hopping is included in the kinetic part 
of the Hamiltonian. The doping level is (a) x = 6%, (b) 
X = 20%, (c) X — 30%. The spectral weight scales as lin- 
ear size of the symbols. The open (filled) circles are re- 
sults of SMA (TMA). The fiUed-blue, fiUed-red and fiUed-pink 
symbols denote results for quasi-electron states, 1-band from 
states |1) and 2-band from states |2), respectively. 



nodal direction and antinodal directions. For instance, 
in Fig.l(a) Z]"(fc = 0) = 0.06 while (fc - kp) = 0.16. 
The band width below the Fermi level is roughly OAeV. 
In the mean time, we also get a band at higher bind- 
ing energy (1.2 ^ 1.5 eV) contributed by the states |2). 
The suppression of low energy spectral weight near the 
zone center and the appearance of a higher energy band 
are exactly what ARPES reveals. Moreover, the energy 
scale of the higher energy band is of order 3t rather than 
the Mott Hubbard gap (which would be infinite in a t-J 
model), it therefore reflects the information of soft lower 
Hubbard band. [13] 

By inspecting Figs.l we see that with increasing dop- 
ing the spectral weight near (below or above) the Fermi 
level increases, while it decreases on the high energy band 
below the Fermi level. This signifies a gradual transfer 
of the spectral weight from the high energy band to the 
low energy bands. On the other hand, the slope of the 
dispersion near the nodal Fermi point barely changes for 
the 1-band under our TMA. The ratio among the Fermi 
velocities we estimated is 1 : 1.05 : 1.25 for x = 6%, 20% 
and 30%. This is rather consistent with the "univer- 
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FIG. 2: (Color online) The same plot as Fig. 1(a) except that 
a third nearest neighbor hopping is added in the Hamiltonian. 



sal Fermi velocity" observed by ARPES. (Due to poor 
screening in cuprates, long-range Coulomb interaction 
at lower doping levels may have caused a slightly larger 
Fermi velocity on top of the correlation effects under dis- 
cussion.) In contrast, the SMA (open-green symbols) 
gives a corresponding ratio 1 : 1.26 : 1.52, [25] and would 
change more significantly in the renormalized mean field 
theory, [30, 31] which is rather far from the experimental 
result. We notice that the universal Fermi velocity was 
also achieved in the literature by tuning J/t [25] or by in- 
cluding a three-site hopping term in the Hamiltonian. [20] 
Finally we observe from Figs.l that in the range of dop- 
ing levels we studied, the energy scales for the low and 
high energy bands below the Fermi level does not change 
significantly. This is in agreement with the experimental 
results. [3, 5] 

It turns out that the suppression of the low energy 
spectral weight near the zone center can be made more 
complete by including the third nearest neighbor hopping 
integral in the Hamiltonian and the trial wave function. 
For doping level x = 6%, h = OAeV, ts = 0.06ey and 
J = 0.12eV the optimal variational parameter in the trial 
wave function is t^^^ — 0.12, A^ar = 0.55, fi^ar — —0.28. 
The quasiparticle spectra are plot in Fig. 2 (a). In this case 
Zf (A: = 0) = 0.045 and Z^{k ~ kp) = 0.18. The curva- 
ture of the high energy dispersion near the zone center 
is changed as compared to the case in Fig. 1(a). Similar 
tendency is found in Fig. 2(b) as compared to Fig. 1(b). 
We observe that the curvature also varies by inspecting 
published data for different families of cuprates. [2, 3, 5] 
We also notice that a negative second nearest neighbor 
hopping integral t2 can lead to a less prominent suppres- 
sion of the spectral weight near the zone center, but the 
existence of a higher energy band is robust. 

We have also investigated the quasiparticle excitations 
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FIG. 3: (Color online) Energy dispersion and spectral weight 
map in the energy-momentum space (along high symme- 
try cuts) in an electron doped case with x = 16%. Open 
(filled) symbols show resuhs of SMA (TMA). The fiUed-blue, 
fiUed-red and fiUed-pink symbols denote results for quasi-hole 
states, 1-band from states |1) and 2-band from states |2) (in 
the electron picture), respectively. Notice that the high en- 
ergy anomaly now appears above the Fermi level. 

in electron doped cuprates, which is still described by 
t — J models after a particle-hole transformation. We 
choose a doping level x = 16%, at which it is believed 
that antiferromagnetic order does not exist. The parame- 
ters in the Hamiltonian are ti = —OAeV and J = 0.12eV. 
The variational parameter are the nearest neighbor sin- 
glet gap function A^^^ with d-wave pairing symmetry and 
the chemical potential jivar- The optimal parameters we 
found are A^^^ = 0.27 and fiyar = —0.3. In order to 
relate to the results in ARPES experiments one need to 
switch back to the electron picture. Our result is shown 
in Fig. 3. The filled circles are the results under TMA, 
and the open circles are from SMA. We find that in the 
electron-doped case the quasiparticle excitations below 
the Fermi level exhibit neither suppression of spectral 
weight at the zone center nor a higher energy band. This 
is consistent with the results of exact diagonalization.[14] 
Instead, a high energy band does appear but above the 
Fermi level. In fact, apart from some change of param- 
eters in the Hamiltonian, the electron-doped dispersion 
and the associated spectral weight can be regarded as 
those of the hole-doped case but viewed up-side-down. 

To conclude, we proposed a two-mode approach to en- 
force the sum rules for quasiparticle excitations both be- 
low and above the Fermi level in our Gutzwiller projec- 
tion variational Monte Carlo study of doped Mott insu- 
lators described by the one-band t-J model. The TMA 
resolves a long standing issue regarding the sum rule in 
VMC.[27] In the hole doped case, we obtained results 
in agreement with the highly unusual features revealed 



in recent angle-resolved photoemission spectra: 1) an al- 
most doping-independent Fermi velocity, 2) two energy 
scales below the Fermi level, and 3) suppression of the 
spectral weight near the zone center on the low energy 
band below the Fermi level. In the electron doped case 
we predicted that the dispersion below the Fermi level 
does not have high energy anomalies. 
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